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Symmetry analysis of the Dzyaloshinskii-Moriya (DM) interaction in MnSi-type cubic crystals 
demonstrates that the magnetic moments are tilted periodically, producing a weak antiferromagnetic 
pattern, when the helix is unwound by magnetic field. The tilt angles of four Mn sublattices are 
determined by a component of the DM vector perpendicular to that one responsible for helical 
spiraling; both components have been evaluated using a simple model. It is shown that the tilting 
should induce pure magnetic reflections 00£, I — 2n + 1 in neutron or x-ray magnetic scattering, and 
the structure factors of these "forbidden" reflections are calculated for arbitrary field orientations. 

PACS numbers: 75.25.-j, 75.10.Hk, 61.05.fm, 61.05.cp 

The crystal magnetism is closely related to crystal structure but the details of this relation are not obvious even for 
simple cases. Meanwhile a more clear understanding of microscopic interactions behind numerous types of magnetic 
ordering would be crucial for magnetic crystal engeneering and applications. Interplay between crystal structure 
and magnetism becomes even more intriguing in noncentrosymmetric crystals where two other players, chirality and 
ferroelectricity, can be also very important. A bright example is the Dzyaloshinskii-Moriya (DM) interaction [J-Q: 
being simple in form, it can result in very different magnetic phenomena: weak ferromagnetism in rhombohedral oxides 
[f| , multiferroics and exotic magnetics in rare-earth compounds @, @] , helical and Skyrmion structures in MnSi-type 
crystals etc. It is well understood that this interaction is related with asymmetry of atomic arrangements. 

However, in spite of numerous theoretical and experimental studies, we have only qualitative understanding of these 
phenomena rather than quantitative descriptions from the first principles. To move further it would be important to 
develope new quantitative methods for experimenal studies revealing the physics of underlying magnetic interactions. 
Recently we suggested [ljj and realized [12( a method for experimental determination of the sign of the DM vector 
in the weak ferromagnetics, and compared the results for FCBO3 with ah initio simulations. 

In this paper we show that the picture of magnetic ordering in non-centrosymmetric crystals of the MnSi type is 
more rich than it was thought before. Even in a strong magnetic field when the helical order disappears, there is still 
a periodic tilting of magnetic moments induced by some component of the DM vector. We estimate the value of this 
component from a simple model and discuss how it can be measured using x-ray or neutron diffraction. 

To find a magnetic pattern formed by four spins s 4 in a unit cell, the classical-spin model, frequently accepted for 



ur spins s 

the magnetic energy E with the DM term [5|-|7|, LL3( , will be used in the following form 

E = E (- Js * • sJ + D " • [ gt x sJ 'D - ^ H • st }> (!) 

where J > is the exchange interaction constant, D* 3 are the DM vectors of the bonds connecting each of four spins 
s* with its six neighbors, and fj, s is an effective magnetic moment of Mn atoms. 

For zero magnetic field the ground state of the system is a spin helix which can be reorientated rather easily and 
becomes conical in not very strong fields. When the helix is unwound by a strong enough external magnetic field 
(|H| > H c2 @, EH): the moments are aligned periodically but, contrary to conventional ferromagnets without DM 
interaction, the moment directions s* may be slightly different for four different types of equivalent atomic positions 
in the MnSi unit cell [IBj . For simplicity the temperature dependence of the moments is not taken into account, all 
the moments are normalized so that |s*| = 1. The upper index t — 1,2,3,4 enumerates four atomic positions in the 
(x, x, x), r 2 = (—2;, I + x, I — x), r 3 = (-| — x, —x, \ + x), r 4 = (^ + x, i — x, —x), correspondingly. 
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Supposing that the DM coefficient is smaller than exchange one (ID* 3 | J), we can assume that s* are tilted from 
s = H/|H| by small angles <p* _L s, and express the energy of the unit cell of the crystal as 



(2) 



{ij} 



i=l 



where the first summation {ij} is taken over twelve differently directed bonds with DM vectors D™ listed in Table 
U and (f l; > = tpi — ip l . 

TABLE I: Twelve bonds between neighboring Mn atoms and corresponding Dzyaloshinskii-Moriya vectors. The bond ij is 
directed from the atom of type i to the atom of type j. An arbitrary DM vector is assumed for bond 12 and all other bonds 
(and their DM vectors) are obtained by symmetry transformations of the cubic point group 23. In MnSi crystal x = 0.138. 
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(-2x, 
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(D x , —D y , —D z ) 



The minimization of ^ gives 



4J + fi s H/2 



(3) 



where vectors are directed along threefold axes passing through r* and interrelated by the symmetry transformations 
of the space group P2i3: 



tA 1 = (1,1,1), 

V 2 = (-1,1,-1), 
V> 3 = (-1,-1,1), 

^ 4 = (1,-1,-1). 



Thus, in the first approximation, the four spins are: 



s* = s + x s], 8 = 



D x + D z 
4J 



(4a) 
(4b) 
(4c) 
(4d) 



(5) 



where the magnetic term is neglected in comparison with exchange, fj, s H <C 8 J. From this equation one can see that 
four components [if)* x s] are always aligned antiferromagnetically, so that X)t=i[' i /'* xs ] =0 f° r an y field direction. 
They give no contribution to the macroscopic magnetization, but they can be measured by diffraction. In particular, 
they induce pure magnetic reflections in x-ray or neutron magnetic scattering. 

Indeed, let us consider neutron scattering (elastic and coherent). It is determined by the structure factor which is 
a sum of nuclear and magnetic structure factors F(h) = F nuc (h) + F mag (h) where h = 2n(h, k,£) is the scattering 
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vector (in our case it is the reciprocal lattice vector with the Miller indices hkt). For simplicity it is supposed that 
magnetic moments are centered on Mn atoms and possible itinerant nature of MnSi magnetism is ignored. In this 
case F mag (h) can be written as a sum over all local moments with corresponding phase multipliers 



4 

F mag (h) = b m P • V (s* - h(h • s^/h 2 ) exp(ih • r ( ) (6) 



t=l 



where P is the direction of neutron spin, b m is the amplitude of neutron magnetic scattering which is proportional to 
/i s , to magnetic formfactor / m (h) and to Debye- Waller factor Jdw- Taking into account that in the magnetic field 
P is either parallel or antiparallel to the field direction so that P = P • s = ±1 and P • [ip* x s] = we obtain 

F m ag (h) = FpM + FafM 
4 

= b m P (1 - (h • s) 2 /h 2 ) ex P(* h ' r ') ( 7a ) 

i=l 

h s 4 

+Sb m P-^ r [h x s] • ^i/>*exp(ih • r*). (7b) 
t=i 

We see that the antiferromagnetic contribution to the structure factor, Eq. (0}), induced by the DM interaction, 
vanishes if the direction of magnetic field is either parallel or perpendicular to the scattering vector whereas the 
ferromagnetic contribution, Eq. ((7^), as it is well known, vanishes if the field is parallel to h. However, owing to 2i 
screw axes, both ferromagnetic and nuclear contributions to the structure factor should disappear for special type 
of reflections, 00£,£ = 2n + 1 (forbidden reflections) because for them J2t=i ex P(^ ' r *) = 0- On the contrary, the 
antiferromagnetic contribution (JTJd) , proportional to X)t=i V'* exp(ih-r'), does not vanish for these reflections and the 
structure factor is given by 

F(00£, £ = 2n+l) = F AFM (00£, £ = 2n + I) 

= Ab m P5sz(—Sy cos 2it£x + is x sin 2n£x). (8) 

The intensity of these reflections, which is proportional to ^(h)! 2 , can be strongly changed by rotation of the field. 
It is easy to prove that for allowed reflections of this type, 00£, £ = 2n, the antiferromagnetic part gives no contribution 
because vector Ylt=i ^ ex P(*h • r *) is parallel to h. For oblique directions of h and s all the contributions to the 
structure factor are nonzero and interference between them allows us to measure the absolute value and the sign of 5. 

There is another phenomenon which can contribute to the forbidden reflections in MnSi, the anisotropy of atomic 
susceptibility tensors Il7j |. previously observed in Nd 3 _ x S4, Sm 3 Te 4 , and i? 2 Ti 2 07 crystals (i?=Ho, Tb, Er, Yb). 
The magnetic susceptibility of each atom can be considered as a second-rank tensor Xjk obeying the point symmetry 
of the atomic site. The magnetic field induces the atomic magnetic moment Mj = XjkH^ even in the paramagnetic 
phase. For rare earths the values of the induced moments can be of about fis in fields of about one Tesla [UBS; for 
transition metals they are never measured. 

It is important that in complex structures the induced moments have different orientations for different atomic sites 
and can contribute to the neutron scattering amplitude of forbidden reflections. In MnSi, where Mn atoms are at 
sites with threefold symmetry axes of four different orientations if} , the site susceptibility tensor can be written as 

X t jk =XoSjk+Xa(i> t M- s jk) (9) 

where 5jk is the unit tensor, xo is the isotropic susceptibility whereas Xa describes the anisotropy. The induced 
magnetic moment is (xo ~ Xa)H + Xa^iip* ■ H) and after simple calculations we obtain the following contribution to 
the structure factor 

F xAFM (00£,£ = 2n + l) 
= 4:b m Ps z (s y cos 2tt£x + is x sin 2ir£x). (10) 

where b m is now proportional to XaH instead of jj, s . Thus, contrary to Fafm, the absolute value of F x afm is growing 
proportional to H, and this allows us to distinguish these two contributions to the reflection intensity. 



X-ray scattering is also sensitive to spin distribution 18-20] and can be applied instead of neutrons to this type of 



studies. It has some important advantages, for instance the polarizations of incident and diffracted beams, e and e', 
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are changeable regardless the direction of the magnetic field, and, in addition, interference with resonant scattering 



21j can be quite helpful. In our case of forbidden reflections from MnSi-type crystals, two contributions to the x-ray 
structure factor are most relevant, nonresonant magnetic scattering and resonant dipole-dipole scattering. The latter 
is expected to be stronger than the former near the if -edge but practically disappears far from the edge. Thus, 
changing the x-ray energy fiw, we can select the best conditions for interference between the two contributions. 

The dipole-dipole tensor of atomic scattering factor looks exactly like Eq. ©: fj k = foSjk + fafyjipk ~ $jk)- F° r 
the "forbidden" reflections OOi', I = 2n + 1 the resonant dipole-dipole scattering changes polarization from tr to it' or 
from 7r to <t [22| and we consider only these diffraction channels 

F a7T , {00e, l = 2n + l)=F dd + Fafm 

4 

= /„ ■ tfW • tr) exp(zh ■ r<) (11a) 

t=i 

h 4 

+iS t i s f m (h) — [s XB]-W exp(ih • r ( ), (lib) 
mc z * — ' 

t=i 

where B = [tt' x tr] + [k' x 7r'](k' • tr) — [k X <r](k ■ n') - [k' x n'] x [k x tr] = 2k sin 2 6, 6 is the Bragg angle, k and 
k' are the unit wave vectors of incident and diffracted waves. We consider only spin scattering [23j | whereas orbital 
scattering is neglected. Finaly we obtain 

F a7I , (00£, £=2n+l) = F dd + F AFM 
= 2f a sin 2#(cos (p cos 2ir£x — i sin ip sin 2ir£x) (12a) 

— 2<5/i s / m (h) sin 2 9[(s x sm6 + s z cos#cos</?) sin27rfo 

mc 2 

+i(s y sin 9 + s z cos 6 sin ip) cos 2ir£x] , (12b) 

where tp is the azimuthal angle of k, starting from the ir-axis. We see that both F dd and Fafm have pronounced 
azimuthal dependence and Fafm can be changed by changing s. There are quite reliable codes for calculation of 
energy dependence of f a , for instance 13, 25|, so that one can, in principle, obtain both the value and the phase 



of S from interference between F dd and Fafm- If magnetic moments have a more complicated distribution, as it is 
expected for itinerant magnetics, Fafm and F x afm could be quite different from those ones given by Eqs. (|T2b ') and 

We consider in detail only forbidden reflections, but their intensity is expected to be rather small. Perhaps it would 
be more appropriate to refine the magnetic structure using many nonforbidden reflections, as it has been done in 
HI 12 1) because from this refinement one can deduce not only the absolute values but also the signs of D x + D Zl 



Xa, and f a . Besides, the forbidden reflections can be excited owing to magnitostriction: for oblique direction of the 
magnetic field the space group reduces from cubic to triclinic, the atoms change slightly their positions, and previously 
"forbidden" reflections become allowed. 

In the absence of ab initio calculations, we can estimate a possible value of the DM vector from a simple expression 
suggested many years ago by Keffer [26|. Starting from the Moriya theory Q of superexchange coupling via a 
near-neighbor intermediate nonmagnetic atom, he had shown that the DM vector could be written as 

D = 0[r 1 'xr 2 '], (13) 

where D is an unknown coefficient. As above, this DM vector corresponds to the Mn-Mn bond directed from r 1 
to r 2 and vectors r 1% and r 2 * are directed from r 1 and r 2 to the position of the intermediate atom. In the MnSi 
structure, there are three intermediate Si atoms per each Mn-Mn bond but two of them form an almost perfect 
parallelogram with two Mn atoms. Therefore their contributions are mutually compensated and we should consider 
only the contribution of one asymmetric Si atom. For this bond it is positioned at (—xsi, \ + xsi, \ — x$i), and Eq. 
(fT3| gives after simple calculations 

D = l -D\x - x Sl \ (2 - 4|a;|,l,4|a:| - 1) , (14) 



This equation is written in the form valid for the MnSi structures of the opposite chirality which are related by 
simultaneous inverting the signs of x and xsi ■ 
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We can find the value of D from the ma gnet ic helix period which is about 40 unit cells (180 A); thus, the helix wave 
vector k « 2tt/40 » 0.157. It was shown [15] that k is related with the DM vector as k = 2\D X - 2D y - D Z \/(3J); 
thus, according to Eq. (fl"4"|). k — \D(x — xsi){l — 8|x|)|/(3J). Accepting x — 0.138 we have \D\\x — x$i\/ J — 4.53 
and \8\ = \D X + £> Z |/(4J) = |-D||x — Xsi\/(8J) « 0.57. We see that parameter 5, which determines the tilt angles, 
is not very small even in MnSi where the helix period is quite long (for comparison, in MnGe it can be about six 
unit cells [27} ). Of course, so large value of 5 seems to be a result of our simple model, Eq. (fT3"|). forcing the 
DM vector to be perpendicular to the bond and ignoring the itinerant nature of magnetism in MnSi. However, 
the DM vector is perpendicular or almost perpendicular to the bond in more sophisticated calculations 0, 0, 0- 
l3lj . Moreover, according to Eq. (fT4|) . the third component of the DM vector, D x + D y — D z is unnaturally large, 
\D X + D y — D z \/J = 2\D\\x — X5i|(l — 2[x|)/J = 6.56. To clarify all those points it would be very important to have 
ab initio calculations of the DM vector for MnSi-type structures. 

In conclusion we have demonstrated that the picture of magnetic ordering in MnSi-type crystals is nontrivial 
even after unwinding of the helix by external field. The remaining antiferromagnetic ordering can be measured by 
diffraction methods providing information about the component of the Dzyaloshinskii-Moriya vector never measured 
before. Our naive estimation shows that this component can be quite large and this is a challenge both to experiment 
and to more reliable calculations. 
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